A Duality Between Unidirectional Charge Density Wave Order and Superconductivity 
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This paper shows the existence of a duality between an unidirectional charge density wave 
order and a superconducting order. This duality predicts the existence of charge density wave 
near a superconducting vortex, and the existence of superconductivity near a charge density wave 
dislocation. 



It is well known that depending on the doping level the 
high T c superconductors exhibit antiferromagnetism and 
d-wave superconductivity. Tranquada et al showed that 
for the high T c compound La2- x Sr x CuO±, four incom- 
mensurate neutron peaks appear around (ir, n) when ap- 
propriate amount of Sr are replaced by Nd . The po- 
sitions of these peaks are consistent with the earlier the- 
oretical proposal of an unidirectional spin density wave, 
or stripe, order 

Recently Lake et al showed that in the highest T c 
LSCO a moderate magnetic field can considerably lower 
the energy of the incommensurate peaks |||. One in- 
terpretation of this result is that the stripe order closely 
competes with superconductivity so that a small suppres- 
sion of the latter by magnetic field allows the former to 
appear (4|J|]. Based on this interpretation a number of 
interesting theoretical works appear recently. In these 
works the LSCO system is viewed as having two compet- 
ing order parameters - the stripe order and superconduc- 
tivity @. 

The purpose of the present work is to extend the 
Landau-like theories in Ref. [^|to include quantum fluc- 
tuations of the stripe and superconducting order param- 
eters || . This theory is relevant when the Landau theory 
predicts the coexistence of both orders. A caveat of the 
present theory is that it assumes an unidirectional charge 
density wave instead of spin density wave order. For this 
reason this theory is not obviously applicable to the ex- 
periments reported in Ref. g]. 

Although this theory is motivated by the "stripe phe- 
nomena" in high T c , it is should be perceived more gener- 
ally as a quantum theory for competing superconducting 
and unidirectional charge density wave orders. Since this 
theory is phcnomcnological in nature it sheds no light on 
the superconducting pairing mechanism and the origin 
of stripe order. Our main results are summarized in the 
last paragraph of this paper. 

The basic assumption we make in this work is the 
coexistence of the stripe and superconducting orders in 
the mean-field theory. Our goal is to go beyond mean- 
field and describe the soft collective fluctuations. Unless 
stated otherwise the entire discussion assumes zero tem- 
perature. The effective actions presented in this work are 
the low-energy and long wavelength quantum Boltzmann 
weight in 2 space and 1 (imaginary) time dimensions. 
The space cutoff of these actions is the maximum of the 



superconducting paring length and the distance between 
neighboring stripes in the mean-field theory. 

We begin with the theory for superconducting phase 
fluctuations 0. When the modulus of the supercon- 
ducting order parameter ip sc is well-developed, the low- 
energy and long-wavelength excitations are associated 
with the fluctuations of the superconducting phase 9 S 
(ipse — \ipsc\e 19 ")- Under that condition the following 
U(l) sigma-model describes the physics at long distance 
and low energy 
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In the above <j> — e l9s is a U(l) phase factor. The period- 
icity of 4> under 9 S — > 9 S + 2mr is crucial for the existence 
of vortices. The last term of the above equation is com- 
monly referred to as the Berry's phase j7j. R$ acts as 
a chemical potential that drives the Cooper pair density 
away from integral number per unit cell. By introducing 
auxiliary fields J x ,J y ,Jt we can Hubbard-Stratonavich 
transform Eq. ([H) and rewrite it as 
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-Jl + J a 



Rd 



(2) 



We note that after integrating out J M and using the iden- 
tity ((fidftcf)) 2 = -|9 M 0| 2 Eq. (0) is recovered. 

If we separate 9 S into a part with winding (i.e vortic- 
ity) and a part without, and integrate out the latter we 
obtain the continuity equation |?J ^2 d^J^ = 0. Phys- 
ically J M is the 3-current of the Cooper pairs. Due to 
the continuity equation we can find a "gauge field" et M so 
that J M = y2u\ e nv\d v a\. Substitute the above equation 
into Eq. (0) we obtain 
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(3) 



In the above f u = X^a e ^\d v a\ is the dual field strength 
of <2a ; aex,/i = (0, a ex ,i, <W,2) with d\a ex a - d%a ex ,\ = R4 
and Va = -i J2„\ ^ux9, 



is the vortex 3-current 



Eq. ([|) describes the physics of Eq. ([!]) from the point of 
view of vortices. The Cooper-pair current J a becomes the 
electric and magnetic fields that couple to the (bosonic) 
vortex current minimally, R^ becomes a background 
magnetic field. 
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Eq. (|l|) involves four parameters iC/>o,i,2 and R^. When 
Rcf, — Eq. ([]]) is the anisotropic 3D xy model. Depend- 
ing on the strength of Ka^ the system can be supercon- 
ducting or insulating. In the superconducting phase the 
vortex current V p are absent at low energy /long wave- 
length. Thus one can view the superconductor as a vor- 
tex insulator. In the insulating phase the vortex cur- 
rent fluctuates severely (i.e. vortices become supercon- 
ducting) which quenches the Cooper current fluctuation. 
Thus the Cooper pair insulator is a vortex superconduc- 
tor. When ^ the vortex condensation is frustrated 
by a background magnetic field. In that case Eq. (|l|) 
describes a superconductor. 

Next we turn to the theory for fluctuating unidirec- 
tional charge density wave. A review of the classical 
statistical mechanics of fluctuating density wave order 
parameters can be found in Ref. [p[. In the following we 
generalize the treatment of Ref. [pf to include the imagi- 
nary time dimension. 

For concreteness let us assume the stripe order causes 
the charge density to modulate at an incommensurate 
wavevector G = Gx. The Fourier component of the 
charge density at this wavevector serves as the order pa- 
rameter. Since we assume the breaking of rotational sym- 
metry the "star" of the ordering wavevectors consists of 
a single pair G and — G. This leads to a single complex 
number ip c d w (and its complex conjugate ip* dw ) as the or- 
der parameter. We further assume the modulus of 4> c dw 
(which is a complex scalar is well developed. The low 
energy excitations are dominated by the fluctuations in 
the phase 9 p of ip cdw (i.e. ^ cdw = \^ c dw\e l9 '')- Physically 
the Op fluctuation is due to the displacement of stripes 
from their ideal position. If u(r, t) is the displacement 
field, the relation between 9 P and u is 9 P = G • u. 

Let us define p = e p , the low energy and long wave- 
length stripe dynamics is described by the following ef- 
fective Lagrangian 

£p= E ^fl^l 2 . (4) 

fi—x,y,t 

The p — x,y terms of the above equation are the elas- 
tic energy due to the stripe displacement, and the p = t 
term is the associated kinetic energy. 

Similar to the passage from Eq. (]]]) to Eq. ^ we 
Hubbard-Stratonavich transform Eq. (^) by introducing 
auxiliary fields q p 

C p = X^or^m + 1»P d M- ( 5 ) 

By separating 9 P into a part involves dislocations and a 
part does not and integrating out the latter we obtain 
^2fi ^nln = 0- Physically qo is the momentum density 
and gi^ are the components of the momentum current 
associated with the stripe displacement. Due to the con- 
servation of q p we can find a "gauge field" b p so that 



Qfi = £nv\dvb\. Substituting the above into Eq. (^) we 
obtain 




In the above f p is the dual field strength of b p and 
D p = —it pv \dvpd\P is the 3-current of the bosonic dis- 
locations. Here a few words should be said about D p . 
In general Dq for a two-dimensional charge density wave 
should be a vector. The fact that Dq in the present the- 
ory is a scalar is due to the one-dimensional nature of the 
charge density wave. 

Depending on the strength of K pp Eq. ([!]) can describe 
the stripe-ordered or the stripe-disordered phases. In the 
stripe-ordered phase the fluctuations of D p are absent at 
low energy and long wavelength, while in the stripe dis- 
ordered phase q p fluctuation is absent. Thus the stripe 
ordered phase is an insulator of the dislocations, while 
the quantum disordered state of stripe can be viewed as 
a Bose condensate of dislocations. 

In the above we assumed that the charge density wave 
ordering wavevector G is incommensurate with the recip- 
rocal lattice vector of the underlying lattice. When G is 
commensurate with the reciprocal lattice vector pinning 
becomes an important issue. The lowest order effect of 
pinning is to exert the following potential on the displace- 
ment field u(r) (recall 9 p (r) = G • u(r)) 

Lpin = ~r E cos ( Gx i + 9p( r i))- ( 7 ) 

i 

Here i label the lattice site and is the corresponding 
position vector. When the stripe period is p/q times the 
lattice constant Gxi — {2-Kq/p)mi where m; is an integer 
between and p—l. Unless p — q = 1 this site-dependent 
phase leads to cancellation when we perform the lattice 
sum. The lowest order pinning term that survives the 
lattice summation is 

L' pm = -r'J2cos( P e p (r l )). (8) 

i 

Strictly speaking such term always pin the stripe at zero 
temperature (in two space dimensions) || . However it is 
important to bear in mind that for large p the pinning 
strength T' ~ T p can be very small, and hence the pin- 
ning is only important at extremely low energies. In that 
case for practical purpose we can ignore the pinning. 

Having understood the effective theories for the su- 
perconducting and the stripe orders separately, we are 
ready to treat them simultaneously. The most impor- 
tant question is how to coupled theses two order param- 
eters together. To answer this question we take hints 
from experiments. Empirically it is known that changing 
the doping density affects the stripe period. In partic- 
ular for sufficiently low doping the stripe period is in- 
versely proportional to the doping density 113]. This 
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suggests that a non-zero Jt tends to shift the phase of 
p, which in turn means a change of 



p, i.e. p 



igiJ t x 



the stripe ordering wavevector from Gx to (G + g\Jt)x. 
Similarly a non-zero current J x tends to boost the stripes 



and shift the phase of p according to p 



p. This 



shift changes the stripe displacement in a time-dependent 
fashion: u(r, t) — ► u + g2J x tx- These considerations sug- 
gest that Eq. (g) and Eq. ([|) be coupled together as follow 



■\{d x - igiJ t )pY 



\{d t - ig2J x )p\ 
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PV \dyP\ 2 



Y}^ J l + + ^ w- ( 9 ) 



We note that in Eq. ([|) Jt and J x act like the x and t 
components of a "gauge field" coupling minimally to the 
gradient of the stripe order parameter, suggesting that 
stripe order tends to suppress the fluctuation in J x j. 

Analogously to the passage from Eq. (||) to Eq. (||) 
we introduce auxiliary fields to Hubbard-Stratonavich 
transform the first three terms of Eq. (0). The result is 



+q x (pdxP - igiJt) + q y pd y p + qt{pd t p - igiJx)- (10) 
Integrate out Jp in the above equation we obtain 



K, 



fm 2 



IK ^ 



q^pd^p] + R^4>d t 4>- (11) 



Comparing Eq. (11) with Eq. (||) we see that aside from 
the last term these two equations can be made identical 
after the following replacements 



p^-> 0, K, 



K„ 



5i ^ 92- 



(12) 



Eq. (|T^) is the main result of this paper, it implies that 
when i?^ = the stripe order and the superconducting 
order are dual to each other. 

Eq. (j^) generically supports four different phases : a) 

< 4> >^ 0,< p >^ ; b) < 4> 0,< p >= 0; c) 

< 4> >= 0, < p >^ ; and d) < 4> >= 0, < p >= 0. 
Since both disordered phase is not terribly interesting, in 
the rest of the paper we shall concentrate on a)-c). 

Phase a) has both <fi and p ordered. (It is important 
to note that by < <j) >^ 0, < p >^ we refer to the co- 
existence of the superconducting and stripe order above 
the cutoff length, i.e., the stripe period. It does not im- 
ply the homogeneous coexistence of these two types of 
order at microscopic scale.) In this phase we can ig- 
nore the vortices in <f> and the dislocation in p and arrive 
at a quadratic theory with two linearly dispersive gap- 
less normal modes. It is straightforward to compute the 



density-density and current-current correlation functions 
using this quadratic theory. The results are summarized 
as follows. 



Kd, x cos 2 9 + Ksy sin 2 i 



K$ x cos 2 9 + (K^y + g(K<f, x K<foKpt) sin 9 



7T q = K^t 



K px cos 2 9 + K py sin 2 9 

K px cos 2 9 + (K py + glKpxKpyK^t) sin 2 9 



Here k w is the transverse current-current correlation 
function in the limit q — > and u> — > while |q|/cj — > 0. 
Physically it gives the superfluid density measured in 
optical conductivity. n q is the current-current correla- 
tion function in the limit q — * and lu — > while 
uj/\q\ — > 0. Physically it gives the the superfluid den- 
sity in the DC magnetization measurement. Finally 7r 9 
is the DC density-density correlation function. 

Due to the anisotropy the response function depends 
on the angle 9 between the vanishing wavevector and 
x (the stripe ordering direction). We note that in the 
presence of stripes n u ^ K q . In principle this inequality 
can be seen by comparing the optical conductivity result 
with the finite wavelength DC magnetization measure- 
ment. Despite this slight complication the above results 
imply that phase a) is an anisotropic superconductor. 

Next we focus on phase b) where there is superconduc- 
tivity but no stripe order. In this case this theory has an 
interesting implication on the p order in the neighbor- 
hood of a superconducting vortex. Around a supercon- 
ducting vortex the phase of <f> winds by 2n. In order to 
minimize the gradient energy in the first term of Eq. (|ll| ) 
a non-zero qt is induced: 



< <?t(x,i) >= i? p (x) =-i< <j)d x <j) > /g 2 



(14) 



Through the q^pd^p term in Eq. ( |ll| ) this gives rise to 
Rp(x)pdtp. Aside from the spatial dependence of R p this 
looks exactly like the Berry phase term in Eq. (nh. If we 
view p as the phase factor of bosons and Eq. (0) as the 
Lagrangian for a system with integral number of bosons 
per unit cell, then the p disordered phase corresponds 
to the bosonic Mott insulator. In this correspondence 
Eq. ([IJ) acts like a spatial-dependent chemical potential. 
In regions where |i? p (x)| exceeds the Mott gap an ex- 
tra boson density will be induced on the Mott insulating 
background. As the result a boson condensate puddles is 
nucleated. Far away from the vortex |i? p (x)| falls below 
the Mott gap consequently the Mott vacuum remains in- 
tact. This means that stripes exist in the vicinity of a 
superconducting vortex. We note that in the above mech- 
anism it is the screening current of the superconductor 
that induces stripe order. This aspect is similar to Ref. 
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Finally we look at phase c) where there is stripe order 
without superconductivity. The lack of 4> order requires 
R$ — 0. Thus the boson density satisfies the require- 
ment that there is an integer number of Cooper pairs per 
stripe unit cell. In this phase if a dislocation is nucleated 
somewhere, the phase of p will wind by 27r around it. 
Through the same mechanism discussed in the previous 
paragraph a spatially varying < Jt(x, t) >— i?^(x) will 
be induced, which leads to superconducting puddles near 
the dislocation. 

In the presence of disorder we expect the existence 
of a finite density of dislocations. The superconducting 
patches associated with these dislocations form a random 
granular array. It is possible that with Josephson tun- 
neling superconductivity can be sustained among these 
grains. In this connection we noted that in a quite dif- 
ferent context Zaanen and Nussinov made an interesting 
conjecture that the high T c superconducting state corre- 
sponds to the condensation of topological excitations of 
an ordered stripe state jj^ ]. 

In the above we have assumed that the charge density 
wave order spontaneously breaks the 90 degree rotation 
symmetry of the underlying lattice. It is interesting that 
removing this assumption has no impact on the duality at 
all. After adding the orientation order parameter Eq. (|^) 
becomes 



K, 



P-L 



\(d ± -i gi j t ) P f 



K, 



pi 



K. 



+ 
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l 



2K, 



2K, 



2K, 



-J 



rf>t 



+J\\4>d\\4> + J±4>d_i_4> + J t (j>d t 4> + R4,4>d t (j) 



+— \dxx\ 



— \dyXl 



^\d t x\ 2 + V(9). 



(15) 

cos 9d v — sin 9d x where 9 



Here — cos9d x +sm9d y , d\\ 
is the angle between the stripe ordering direction and the 
crystalline axis, x — el9 1S the orientation order parame- 
ter, and V(9) = —K cos 49 acts to pin 9 to 0, 7r/2, tt, 3tt/2. 
Starting with Eq. (||) we can redo duality transforma- 
tion between p and <fi with \ unaffected. 

Doped Mott insulator often exhibit inhomogeneous 
structure of competing orders. Unidirectional charge 
density wave and superconductivity seems to compete 
in the LSCO family of the cuprates . Similar com- 
petition between a charge ordered state and ferromag- 
netism exists in the colossal magneto-resistive materials 
p3[ . In this paper we present an effective field theory 
for competing superconducting and unidirectional charge 
density wave orders |I3]. We find an interesting duality 
fl5[ between these two types of order (Eqs. | ll|l2]). An 
implication of this duality is the presence of stripes near 
a superconducting vortex, and the presence of supercon- 
ductivity near a stripe dislocation. The evidence of the 



former is already obtained in Ref. |3j. To check the lat- 
ter one can apply stress to stripe ordered state to induce 
dislocations. 
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